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PREFACE  AND  SUMMARY 


This  Note  describes  in  mathematical  terms  the  dynamic  geometry 
between  a  constellation  of  satellites  deployed  for  ballistic  missile 
defense  (BMD)  and  the  missiles  the  satellites  are  engaging.  Formulas, 
readily  translatable  into  computer  code,  are  given  for  such  engagement 
parameters  as  slant  range,  closing  velocity,  and  1 ine-of-sight  incidence 
angle  in  terms  of  satellite  and  missile  position  and  velocity  data. 

These  formulas  are  the  foundation  of  a  model  developed  at  The  RAND 
Corporation  to  support  study  of  the  BMD  capability  of  various  satellite 
armament  concepts. 

This  work  was  initiated  at  RAND  under  sponsorship  of  the  Directed 
Energy  Office  of  the  Defense  Advanced  Research  Projects  Agency. 
Sponsorship  was  transferred  to  the  Directed  Energy  Office  of  the 
Strategic  Defense  Initiative  Organization  when  that  agency  was  formed. 
The  study  was  carried  out  in  RAND's  National  Defense  Research  Institute, 
a  Federally  Funded  Research  and  Development  Center  sponsored  by  the 
Office  of  the  Secretary  of  Defense. 
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THE  GEOMETRY  OF  A  SATELLITE-BALLISTIC  MISSILE  ENGAGEMENT 


Our  objective  is  to  describe  the  dynamic  geometry  of  an  engagement 
in  which  a  constellation  of  satellites  is  defending  against  ground-based 
ballistic  missile  launches.  The  formulas  presented  allow  one  to  keep 
track  simultaneously  of  satellite  motion,  missile  motion,  and  earth  rota¬ 
tion,  and  thus  permit  the  determination  of  such  variables  as  slant  range, 
closing  velocity,  line-of-sight  visibility,  and  missile/satellite  orien¬ 
tation. 

ASSUMPTIONS 

All  satellites  are  assumed  to  be  in  circular  orbits  at  a  common  al¬ 
titude.  We  assume  a  spherical  earth  of  uniform  density  and  ignore  possi¬ 
ble  orbit  perturbations  due  to  extraterrestrial  bodies.  The  constella¬ 
tion  consists  of  m  orbit  rings  with  n  satellites  per  ring  and  is  called 
an  w  X  n  aonstellation.  All  rings  have  the  same  inclination  with  respect 
to  the  equatorial  plane,  with  satellites  in  a  given  ring  equally  spaced 
and  traveling  in  the  same  direction.  We  consider  two  ring-spacing  op¬ 
tions  : 

Option  A:  The  ascending  nodes  of  the  rings*  are  equally  spaced 
around  the  equator.  In  the  case  of  polar  orbits  with  m  even,  each 
ring  is  coincident  with  another  ring,  with  satellites  in  the  two 
rings  traveling  in  opposite  directions. 

Option  B:  The  ascending  nodes  of  the  rings  are  equally  spaced 
around  one-half  the  equator. 


;SI 


k'd 


These  options  are  illustrated  for  two  orbit  inclinations  in  Fig.  1. 

We  select  a  fixed  reference  ring,  called  ring  1,  and  label  the 
remaining  rings  2,  3,  4,  ...,  according  to  the  position  of  their  ascend¬ 
ing  nodes  westward  from  ring  1.  This  reference  ring  can  be  any  of  the 


-  .V  . 


That  is,  the  point  at  which  satellites  in  the  ring  cross  the  equa¬ 
torial  plane  and  ascend  into  the  northern  hemisphere. 
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(VIEWED  FROM  INFINITELY  FAR  ABOVE  THE  NORTH  POLE) 


90®  inclination 


Ring-spacing 
option  B 


60®  inclination 


m=3 


m'4 


Fig.  1  —  Examples  of  constellation  configurations 


m  rings  under  option  A;  under  option  B,  it  is  that  ring  whose  ascending 
node  is  the  easternmost  along  that  half  of  the  equator  containing  all  m 
ascending  nodes.  Select  some  fixed  satellite  in  ring  1  and  agree  that 
it  occupies  position  1  in  this  ring.  This  satellite,  labeled  1-1,  will 
henceforth  be  called  the  vefevence  sate'Ll'Lte.  Label  the  positions  of 
other  satellites  in  ring  1  by  2,  3,  4,  ...,  moving  forward  (i.e.,  in  the 
direction  of  satellite  motion)  from  the  reference  satellite.  In  general, 
satell''*-'?  i-j  occupies  position  j  in  ring  i.  The  positions  of  satellites 
in  rings  2  through  m,  relative  to  those  in  ring  1,  are  determined  by  the 
oonstellat'Lon  phasing  fraction  p  (0  <  p  <  1).  In  particular,  we  assume 
that  satellite  2-1  is  Q  =  360p/n  deg  ahead  (in  argument  of  latitude)  of 
satellite  1-1,  satellite  3-1  is  Q  deg  ahead  of  2-1,  and  so  on.  For 
ring-spacing  option  A,  in  order  that  satellites  in  rings  m  and  1  occupy 
the  same  relative  positions  as  those  in  other  pairs  of  adjacent  rings, 
p  is  restricted  to  the  values  0,  1/m,  2/ra,  ...,  (m  -  l)/m.  This  require¬ 
ment  cannot  be  met  under  option  B  since  satellites  in  rings  m  and  1 
travel  in  opposite  directions;  thus  p  can  assume  any  value  between  0 
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and  1.  So  that  we  can  assign  times  to  various  events,  we  assume  that 
at  some  0-hour  time,  the  reference  satellite  crosses  the  equator  at 
E  deg  longitude  heading  north. 

To  facilitate  our  presentation,  we  introduce  an  earth-referenced 
(rotating)  xyz  coordinate  system,  with  the  origin  (0,0,0)  at  the  earth's 
center.  The  equator  lies  in  the  xy-plane  with  the  positive  x-axis  pas¬ 
sing  through  the  Greenwich  prime  meridian,  the  positive  y-axis  lying 
90  deg  east,  and  the  positive  z-axis  passing  through  the  north  pole. 

As  is  conventional,  latitudes  are  measured  northward  from  the  equator, 
longitudes  eastward  from  the  positive  x-axis,  azimuths  clockwise  from  due 
north,  and  path  elevation  angles  upward  from  horizontal. 

The  ballistic  missile  (which  we  refer  to  as  the  ICBM)  is  sited  at 
a  particular  latitude  and  longitude  and  is  launched  with  a  known  azi¬ 
muth  relative  to  the  above  coordinate  system.  It  is  assumed  that  suf¬ 
ficient  ICBM  trajectory  data  are  available  (e.g.,  from  a  powered  flight 
simulator)  to  establish  ICBM  position  and  velocity  vectors  at  all  times 
of  interest.  An  example  of  such  data  is  shown  in  Table  1  (p.  9)  for  the 
boost  phase  period  of  a  hypothetical  ICBM.  Although,  in  general,  six 
trajectory  elements  are  needed  to  establish  position  and  velocity,  the 
actual  six  used  can  vary  depending  on  the  output  of  the  model  used  to 
simulate  the  ICBM  flight.  Rand's  model  COMET  (see  R-32^0-ARPA,  forth¬ 
coming)  outputs  position  and  velocity  directly.  To  ac<'ommodate  models 
that  do  not  directly  output  these  data,  we  describe  in  the  pages  that 
follow  how  to  calculate  tlu'se  vectors  and  (>ther  relevant  parameters, 
assuming  as  basic  inputs  the  Ih'l  lowing  six  ICBM  trajecto'-v  elements: 


•  ICBM  altitude  •ICBM  speed 

•  ICBM  latitude  •ICBM  path  azimuth 

•  ICB.M  longitude  •  ICBM  path  elevation  angle 

The  formulas  incorporate  the  following  definitions  and  conventions: 

Def inltlons 

So.follite  orlit  (f)  —  the  angle  between  the  equatorial 

plane  and  the  plane  of  the  satellite's  orbit. 
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0°  <  e  <  90°  if  satellite  travels  eastward  in  its  orbit 


e  =  90 


if  satellite  travels  due  north  and  south 


90°  <  0  <  180°  if  satellite  travels  westward  in  its  orbit 


I  CBM  "path  elevation  angle  (5) — the  angle  between  the  ICBM’s  velocity 
vector  and  the  local  horizontal  plane. 

ICBM  path,  azimuth  (X) — the  angle  in  the  local  horizontal  plane  be¬ 
tween  the  projection  of  the  ICBM's  velocity  vector  and  the  vector  pointing 
due  north. 

Satellite  argument  of  latitude  (\{)) — the  (non-negative)  central  angle 
in  a  satellite's  orbit  plane  between  the  ascending  node  and  the  satellite. 


Conventions 


1. 

2. 

3. 


Altitudes  are  measured  from  the  surface  of  the  earth. 


-90°  <  sin“^(  )  <  90°; 


0°  <  cos“^(  )  <  180°. 


-90°  <  tan  ^(  )  < 


90 


Arg(x,y)  is  the  angle  (measured  counterclockwise)  between  the 
positive  x-axis  and  the  ray  joining  (0,0)  to  (x,y).  In  parti¬ 
cular. 


arg(x,y)  = 


tan  ^(y/x) 

tan  ^(y/x)  +  180^ 

90° 

-90° 

undefined 


if  X  >  0 
if  X  <  0 
if  X  =  0,  y  >  0 
ifx=0,  y<0 
if  X  =  0,  y  =  0 


4.  The  dot  product  (•)  of  two  vectors  x  =  (x^,X2,X2)  and 

y  =  ^yi>y2’y3^  is  x  •  y  =  X^y^  +  yi^y^  +  x^y^. 

5.  The  cross  product  (x)  of  two  vectors  x  =  (x^,X2,x^)  and 

y  =  (yi’y2’y3)  is  x  x  y  =  (X2y3  -  X3y2,x3y^  -  Xj^y3,x^y2  -  X2y3^). 


6.  The  length  (] 


I  I  I )  of  a  vector  x  =  (x^,x^,x^)  is 

=  /x^  +  x|  +  X3 


7.  Unless  otherwise  indicated,  all  variables  in  the  lists  and  tables 
that  follow  are  referenced  to  the  rotating  xyz  coordinate  system 
described  earlier. 


v->'i 
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INPUTS 

Constants 

Earth  radius  (km) 

3  2 

Earth  gravitation  constant  (km  /sec  ) 
Length  of  sidereal  day  (sec) 

Earth  rotation  rate  (rad/sec) 


Constellation  Variables 

Satellite  orbit  inclination  (deg) 
Satellite  orbit  altitude  (km) 

Number  of  rings  of  satellites 
Number  of  satellites  per  ring 
Ring  number  of  particular  satellite 
Position  number  of  particular  satellite 
Constellation  phasing  fraction 

Ring  spacing  option 

Longitude  of  referer 
0-hour  (deg) 

ICBM  Variables 


(min  after  0-hour) 


Notation 


r 

G 

F 

01 


0 

h 

m 

n 

i 

j 

P 


ICBM  path  elevation  angle  at  time  t  (deg) 
ICBM  path  azimuth  at  time  t  (deg) 


Y 

6 

X 


Range 
6375.58 
398603 
86164.1 
7.29212  X  10 
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0  <  9  <  180 


0  <  p  <  1 


e 

|l  if  option  A 
(2  if  option  B 

*. 

. 

satellite  at 

:v^^: 

E 

•  % 

.XvS 

after  0-hour) 

q 

••  •-  V 
■  ••  < 

ICBM  launch 

-  > 

t 

cr 

A  1 

t  (km) 

t  (deg) 

a 

A 

7-:^ 

t  (deg) 

B 

r  A 

.yj 


'‘V-'.'i 

•  J 

■  -B  ■  I  ^  I  i*>  I  iJt  I  ^  -A .  ^  -*  •  •V*  -A  kV  ■ 


OUTPUTS 


Notation 


Satellite  orbital  period  (sec) 


2-ir 


(r  +  h)- 


Satellite  orbital  speed' 


60  •  360 


deg/min 


km/ sec 


Angular  separation  between  ring  i  and  ring  1 
(deg  east  from  ring  1) 

360(1  -  1) 
me 


Longitude  of  ascending  node  of  ring  i  at 
time  t  (deg) 


$ 


tu 


60  •  180  ^  ^ 
- h  E 

TT 


Argument  of  latitude  of  satellite  i-j  at 
time  t  (deg) 


360(p(i  -  1)  +  1  -  1)  ^ 
n 


Latitude  of  satellite  i-j  at  time  t  (deg) 

.  -1/  .  ,  .  -X 

sin  (sin  V  sin  0) 

I.oneitiidr-  of  satellite  i-i  at  lime  f  fdeiil 


s 


k 


H 


."-V .v.vjvv  ■>  V '>  V 
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Positlon  vector  of  satellite  1-j  at  time  t  (km) 

X(l)  =  (r  +  h)  cos  L  cos  M 
X(2)  =  (r  +  h)  cos  L  sin  M 
X(3)  =  (r  +  h)  sin  L 


Speed  of  satellite  i-j  at  time  t  (km/r ac) 

-  2zw(r  +  h)  cos  9  +  tij^(r  +  h)^  cos^  L 
Path  azimuth  of  satellite  i-j  at  time  t  (deg) 

o 

arg(z  cos  ip  sin  9,  z  cos  9  -  'jj(r  +  h)  cos^  L) 

Velocity  vector  of  satellite  i-j  at  time  t  (km/sec) 

u(l)  =  c(-  'OS  b  sin  L  cos  M  -  sin  b  sin  M) 

u(2)  =  c(-cos  b  sin  L  sin  xM  +  sin  b  cos  M) 

u(3)  =  c  cos  b  cos  L 

ICBM  position  vector  at  time  t  (km) 

Y(l)  =  (r  +  a)  cos  A  cos  B 

Y(2)  =  (r  +  a)  cos  A  sin  B 

Y(3)  =  (r  +  a)  sin  A 

ICBiM  velocity  vector  at  time  t  (km/sec) 


v(l)  =  y(-COS 

> 

COS 

A 

sin 

A 

COS 

B 

- 

sin 

X  cos 

c 

sin  B 

+  sin  5 

COS 

A  cos  B) 

v(2)  =  y(-cos 

>, 

cos 

6 

sin 

A 

sin 

B 

+ 

sin 

>,  cos 

cos  B 

+  sin  6 

cos 

•A  sin  B) 

v(3)  =  y(  cos 

i 

cos 

<5 

cos 

A 

+  sin 

sin 

A) 

Slant  range 

between 

satellite 

i- 

-j 

and 

ICBM 

at 

t  ime 

t  (km) 

d 

V 


/  (X(l)  -  Y(l))^  +  (X(2)  -  Y(2))^  +  (X(3)  -  Y(3))‘ 


Closing  velocity  between  satellite  i-j  and  ICBM 
at  time  t  (km/sec) 

(u  -  v)  >  (Y  -  X) 
d 


m 

•.•v' 


w 

•  I 

w 


.'a'yv'3 

.'A  VJ 

vVv 
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Angle  between  ICBM  velocity  vector  at  time  t  and 
line-of-sight  vector  between  satellite  i-j  and 
ICBM  (deg) 


-1  /  V  •  (X  -  Y) 


Minimum  altitude  along  line-of-sight  vector  joining 
satellite  i-j  to  ICBM  at  time  t  (km) 


Min  {h,a} 


lx  X  Yl| 
lx  -  Yll 


if  [(X  -  Y)  •  X  [(X  -  Y) 


otherwise 


AN  EXAMPLE  ENGAGEMENT 

We  will  illustrate  the  use  of  these  formulas  by  considering  a  hypo¬ 
thetical  engagement  in  which  a  24-satellite  constellation  is  defending 
against  an  ICBM  launched  from  50°  latitude,  90°  longitude,  and  with  an 
initial  azimuth  heading  of  10  deg.  Trajectory  information  for  such  a 
launch  is  given  in  Table  1.  We  will  assume  that  the  launch  occurs  at 
0-hour,  so  that  the  five-minute  boost  phase  period  shown  in  the  table 
occurs  between  times  t  =  0  and  t  =  5.  Constellation  inputs  are  assumed 
as  follows; 


Satellite  orbit  inclination 
Satellite  orbit  altitude 
Number  of  rings  of  satellites 
Number  of  satellites  per  ring 
Constellation  phasing  fraction 
Ring-spacing  option 

Longitude  of  reference  satellite  at  0-hour 


?  =  75  deg 
h  =  1250.75  km 
m  =  6 
n  =  4 
P  =  2/3 
A  (e  =  1) 

E  =  0  deg 


We  compute 


Satellite  orbital  period 


Satellite  orbital  speed 


s  =  6628.01  sec 
^  k  =  3.26  deg /mil 
I z  =  7.23  km/ sec 


./  V  'J 


POSITION  AND  VELOCITY  DATA  FOR  A  HYPOTHETICAL  I  CBM 
(Boost  Phase) 
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The  positions,  velocities,  and  azimuths  of  the  24  satellites  at  the 
time  of  ICBM  launch  (t  =  0),  shown  in  Table  2,  are  found  by  letting  i  run 
from  1  to  6  and  j  from  1  to  4  in  the  formulas  for  ip  and  ip,  and  then  sub¬ 
stituting  in  the  formulas  for  L,  M,  X,  u,  b,  and  c.  For  the  sake  of  example, 
we  select  satellite  1-2  and  show  in  Table  3  its  positions  and  velocities 
during  the  ICBM  boost  phase.  It  is  constructed  like  Table  2,  except  that 
time  is  varied  between  t  =  0  and  t  =  5,  the  end  of  boost  phase.  Finally, 
in  Table  4,  the  data  of  Tables  1  and  3  are  combined  using  the  formulas 
for  d,  f,  and  D  to  yield  the  crucial  engagement  parameters  shown. 

We  see,  for  example,  that  the  slant  range  between  the  ICBM  and  satel¬ 
lite  1-2  changes  very  little  during  the  boost  phase  (thus  the  relatively 
low  closing  velocity).  This  can  be  seen  geometrically  in  Fig.  2,  which 
shows  the  earth  shadows  of  the  orbit  rings  at  0-hour  and  the  ground  trace 
of  the  ICBM  during  its  boost  phase.  The  fifth  column  of  Table  4  gives  the 
angle  between  the  ICBM  velocity  vector  and  the  line-of-sight  vector  between 
satellite  and  ICBM.  The  absolute  difference  between  this  angle  and  90  deg 
represents  the  inaidenae  angle  (measured  from  the  normal  to  the  ICBM  sur¬ 
face)  that  would  be  encountered  by  a  laser  beam  fired  from  the  satellite 
toward  the  ICBM  and  passing  through  the  ICBM  axis. 

The  last  column  of  Table  4  records  the  minimum  altitude  (above  the 
earth's  surface)  of  points  along  the  line-of-sight  vector  joining  thd 
satellite  to  the  ICBM.  If  this  minimum  is  positive,  the  ICBM  is  visible 
to  the  satellite  above  the  earth's  horizon. 

PERIODICITY  IN  THE  ENGAGEMENT  GEOMETRY 

If  the  altitude  h  of  the  m  n  satellite  constellation  is  chosen 
properly,  the  relative  geometry  between  the  rotating  earth  and  the  con¬ 
stellation  repeats  every  fixed  number  of  days.  This  phenomenon  is  an 
analytic  convenience  in  that  it  bounds  the  time  period  that  must  be 
examined  in  determining  the  effectiveness  of  the  constellation  in  defend¬ 
ing  against  ICBM  launches.  In  the  paragraphs  that  follow,  we  determine 
the  minimum  length  of  such  a  "repetition  period." 

We  assume  here  that  the  ICBM  is  cylindrical  with  its  axis  always 
pointing  along  its  velocity  vector. 
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Table  4 

BOOST  PHASE  ENGAGEMENT  PARAMETERS  FOR 
HYPOTHETICAL  ICBM  AND  SATELLITE  1-2 

(ICBM  launched  at  0-hour) 
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Table  4- 

-cont inued 

Angle 

.Minimum 

T  ime 

T  ime 

between  LOS 

altitude 

after 

after 

Slant 

C los ing 

and  ICBM 

along 

launch 

0-hour 

range 

velocity 

direction 

LOS  vector 

(sec) 

(min) 

(km) 

(km/sec ) 

(deg) 

(km) 

t 

d 

f 
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Fig.  2  —  Example  engagement  geometry— 6  x  4  constellation 


At  the  altitude  h,  each  satellite  makes  c  orbits  per  sidereal  day, 


where 


86164.1 


Imagine  a  point  Q  on  the  earth's  surface  which  at  some  time  0  is  at 
0°  latitude  and  0°  longitude.  Assume  that  the  reference  satellite  (1-1) 
is  directly  overhead  Q  at  this  moment,  and  heading  north.  Suppose  that 
p  days  later,  the  earth  occupies  the  same  position  relative  to  the  con¬ 
stellation  as  it  did  at  time  0.  The  determination  of  the  minimum  pos¬ 
sible  value  for  p  varies  depending  on  the  ring-spacing  option  assumed 
(either  option  A  or  B) . 


Variables  (such  as  s)  that  are  not  defined  here  are  as  defined 
earlier . 


L«*i.  A' A. 


» 


Rins-spaclng  Option  A 

The  point  Q,  p  days  after  time  0,  must  again  be  directly  beneath 
some  north-bound  satellite.  This  forces  p  to  be  a  positive  integral 
multiple  of  1/m,  the  fraction  of  a  day  required  for  the  earth  to  rotate 
from  under  one  ring  to  under  the  one  immediately  to  the  east.  Thus 
p  =  I/m  for  some  positive  integer  I.  The  satellites  in  the  Ith  ring 
to  the  east  of  ring  1  are  shifted  ahead  of  those  in  ring  1  by 
p (m  -  I)/cn  days.  Thus  for  some  satellite  in  this  ring  to  be  at  the 
equator  p  days  after  time  0  requires  that 


J  p (m  -  I) 


for  some  integer  J 


Equating  the  two  expressions  for  p  gives 

_  J _ p(m  -  I) 

m  cn  cn 

Set  p  =  l/m  (recall  that  p  can  only  assume  the  values  0,  1/m,  2/m,  ..., 

(m  -  l)/m).  After  substituting  and  simplifying,  we  obtain 

I(cn  -  £)  =  m(J  -  £) 

To  find  the  mininrum  period  p  (that  is,  the  smallest  possible  value 
for  p),  it  suffices  (since  p  =  l/ra)  to  determine  the  least  positive 
Integer  I  for  which  this  equation  has  an  integer  solution  J.  If  c  is 
not  a  rational  number  (that  is,  the  quotient  of  two  integers),  no  solu¬ 
tion  is  possible.  If  c  is  rational,  express  it  as  the  quotient  x/v  of 
two  Integers,  in  lowest  terms.  It  can  be  shown  that  the  minimal  solution 


gcd(nx  -  £,y,  ray) 

where  the  greatest  common  divisor  of  two  integers  F  and  C  (gcd(F,  (i))  is 
defined  as  the  largest  integer  that  divides  evenly  into  both  F  and 


We  conclude  that  the  minimal  period  p  is  given  by 


^  gcd(nx  -  Jly,  my) 

This  computation  is  illustrated  below  for  four  different  constellation 
configurations.  In  example  1,  the  engagement  geometry  repeats  every  day; 
in  example  2,  every  sixth  of  a  day;  in  example  3,  every  three  days;  and  in 
example  4,  every  fourth  of  a  day. 


Example 

m 

n 

P 

Si 

c 

h 

X 

y 

nx-Hy 

my 

gcd(nx-)ly,my) 

A 

P 

1 

6 

4 

V2 

3 

13 

1250.75 

13 

1 

49  .. 

6- 

1 

1 

2 

6 

4 

Va 

4 

13 

1250.75 

13 

1 

48 

6 

6 

Ve 

3 

3 

8 

0 

0 

14V3 

770.15 

43 

3 

344 

9 

1 

3 

4 

4 

6 

V2 

2 

IAV3 

770.15 

43 

3 

252 

12 

12 

V4 

Ring-spacing  Option  B 

An  m  X  n  constellation  configured  under  this  option  does  not  possess 
the  perfect  syrometry  of  one  configured  under  option  A.  This  lack  of  sym¬ 
metry  results  from  the  existence  of  a  "seam"  between  rings  1  and  m,  with 
satellites  in  these  two  rings  moving  in  opposite  directions,  and  those  in 
other  pairs  of  adjacent  rings  moving  in  the  same  direction.  As  a  conse¬ 
quence,  repetition  periods  of  less  than  one  day  are  not  possible.  If  each 
satellite  makes  c  orbits  per  day,  and  c  =  x/y  is  a  rational  number  in  low¬ 
est  terms  (x  and  y  are  integers)  ,  then  the  mini-mm  repetition  period  is 

A 

p  =  y/gcd(y,  n)  days.  This  formula  is  illustrated  in  the  table  below. 
Note,  in  particular,  that  p  is  independent  of  the  phasing  fraction  p. 
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132/3 

1000.68 

41 
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